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ABSTRACT 

The forced motion of superfluid vortices in shear oscillations of rotating solid neutron 
star matter produces damping of the mode. A simple model of the unpinning and 
repinning processes is described, with numerical calculations of the consequent energy 
decay times. These are of the order of 1 s for typical anomalous X-ray pulsars with 
rotational angular velocity « 1 rad s _1 but become very short for radio pulsars with 
fi w 10 — 100 rad s _1 . The superfluid friction processes considered here may also be 
significant for the damping of r-modes in rapidly rotating neutron stars. 

Key words: stars: neutron-pulsars: general. 
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1 INTRODUCTION 

Although there has been much recent interest in shear modes 
of neutron star oscillation (Blaes et al 1989; Thompson & 
Duncan 1995; Duncan 1998), relatively little attention has 
been paid to the associated damping mechanisms apart from 
coupling with magnetospheric Alfven waves. Blaes et al con- 
sider that stress-induced motion of dislocations may be an 
important process, as is believed to be the case for shear 
waves in the Earth. But nothing is known about the dis- 
tribution of linear and planar defects in the lattice of nu- 
clei at matter densities p below the neutron-drip theshold 
Pnd = 4.3 x 10 11 g cm -3 or about their stress-induced mo- 
bility. At densities in the crust above p n d, it is likely that 
neutron star matter is an amorphous heterogeneous solid, 
analogous with a disordered many-component alloy (Jones 
1999), and again, nothing is known about its likely response 
to periodic stress. Thus there is no basis for even a rough 
estimate of the contribution to shear-mode Q-values from 
this mechanism. 

Blaes et al assume that coupling with magnetospheric 
Alfven waves is the dominant mechanism of damping. These 
authors show that, for a typical frequency of 10 4 Hz, shear 
waves are efficiently reflected from the boundary with the 
magnetosphere and, in the absence of other damping pro- 
cesses, would have a lifetime r s ~ 1 s at a reference mag- 
netic flux density of 10 11 G. The B-dependence of r 3 is given 
in analytical expressions valid asymptotically at frequencies 
small and large compared with 10 4 Hz. These are t s oc B~ 2 
at low freqencies and r s oc _B _4//7 at high frquencies. The 
coupling is therefore an efficient mechanism for energy trans- 
fer from crust to magnetosphere at the high magnetic fields 
assumed in magnetars (Thompson & Duncan 1995) . Sudden 
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movements of the solid crust, caused by Maxwell stresses or 
otherwise, may then lead to magnetospheric excitation and 
possibly observable phenomena. In this context, the present 
paper examines some processes, specific to neutron stars, 
which are more accessible to calculation than stress-induced 
dislocation movement and which might be capable of damp- 
ing shear waves enough to invalidate this assumption. Shear 
waves are reflected many times from the boundary with the 
magnetosphere and possibly from the inner boundary of the 
crust with the liquid core of the star where the matter den- 
sity is p c « 2 x 10 14 g cm -3 . The number of reflections 
must be at least of the order of c s t s /R » 10 2 , where R is 
the neutron star radius and c s ~ 10 s cm s _1 is the shear- 
wave propagation velocity. Reflection from each boundary 
is complicated because it couples shear waves with longi- 
tudinal sound (see, for example, Landau & Lifschitz 1970, 
p. 103) and at p c , may even depend on the dimensionality of 
the solid. 

It is known that longitudinal sound is strongly absorbed 
in the core of a neutron star owing to bulk viscosity pro- 
duced by the strangeness-changing four-baryon weak inter- 
action. The magnitude of the coefficient is very dependent 
on the assumptions which must be made, principally about 
the equation of state at core densities, but the absorption is 
always strong enough, within certain temperature intervals, 
to be significant for shear-wave damping. 

The shear-wave velocity depends on nuclear parame- 
ters according to c s oc ZA~ 2 ^ 3 p 1 J 6 , where Z and A are the 
nuclear charge and mass number, and pjv is the spatially- 
averaged matter density derived from the nuclei alone and 
excluding the superfluid neutron continuum external to 
them. This elementary expression assumes a nuclear effec- 
tive mass equal to the free mass and, in conjunction with 
the nuclear parameters listed by Negele & Vautherin (1973), 
gives a rapid increase of c s with depth at p < p n d but no 
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more than a very slow increase at p > p n£ j. (The latter vari- 
ation seems to be in disagreement with Blaes et al who state 
that c s is a decreasing function of depth at p > p n d-) As a 
consequence of this variation, shear waves will be refracted 
radially outward in all regions, as are longitudinal sound 
waves, but only weakly so at p > p n d- But in a typical 1.4M 
neutron star, the layers of strong refraction at p < p nd have 
depth only of the same order as the wavelength at the rele- 
vant frequency of 10 4 Hz. Even assuming, with Blaes et al, 
that the Maxwell-stress generated events producing shear 
waves occur predominantly within this latter region, it fol- 
lows from this small depth and from the multiplicity of re- 
flections that shear-wave propagation in the neutron drip 
region at p > p n£ j must be considered as a significant source 
of damping. Thus most of the present paper will describe 
an investigation in this context of phenomena involving the 
rapid forced motion of vortices which, for want of any better 
name, we shall refer to as supcrfluid friction. The relevant 
results for longitudinal sound absorption in the liquid core 
are considered in Section 3. 



2 DESCRIPTION OF SUPERFLUID FRICTION 

An elementary description of superfluid friction is given by 
considering a region of solid in which the neutron superfluid 
has rotational angular velocity Q n . The superfluid velocity 
v n is defined (see Sonin 1987) as the average over the unit 
cell of the local vortex lattice. All velocities are referred to 
cylindrical polar coordinates rotating with constant angular 
velocity Q, coincident at some fiducial time with the an- 
gular velocity of the charged components of the star. The 
difference Q n — Q, is always small compared with Q. To the 
extent that steady spin-down of the superfluid occurs, it is 
maintained by radial outward expansion of its vortex lat- 
tice; the Magnus force (m acting on an individual vortex is 
balanced by the force per unit length f generated through 
its motion with velocity — u, perpendicular to its local axis 
and relative to the solid. The quantities f and u have to be 
treated as averages over macroscopic lengths and times, an 
assumption which is satisfactory provided vortex curvature 
on these same scales is negligible. A shear wave of wavevec- 
tor q is represented by a solid matter displacement vector 
£ = expi(q-r — cot), with £q = 0. The velocity relative to 
the local solid is then — u = vl-vi, where the shear- wave 
displacement velocity is v = £ and vi is its component per- 
pendicular to the vortex axis. The velocity of a point on the 
vortex axis is vi. Energy dissipation is obviously changed 
by the displacement velocity, specifically, by the component 
vi. Energy is transferred irreversibly from the shear wave 
itself and, for most wave directions and polarizations, from 
the rotational kinetic energy of the superfluid measured in 
the rotating coordinate system. It is to be expected that 
calculation of shear-wave energy-loss rates depends princi- 
pally on the assumptions made about the dissipative force, 
and so we are obliged to consider these in as much detail as 
possible. 



2.1 The dissipative force 

The assumptions, made here for almost rectilinear vortices, 
are that it is possible to define, averaged over macroscopic 



lengths and times, a force f(tt) per unit length of vortex, 
and that fluctuations from it are not so large as to make its 
use meaningless. The component of this force parallel with u 
is the dissipative force f\\{u). Its properties were considered 
briefly in a previous paper (Jones 2002), but it seems appro- 
priate to give here a more detailed account of what can be 
deduced about it from the interaction of a vortex, comoving 
with superfluid of velocity v„, and a single point defect in 
the solid structure. This leads to the excitation of phonons in 
the solid and, at finite temperatures, of vortex-core quasipar- 
ticles (Jones 1990a, b), but Kelvin wave excitation is easily 
the most important process (Epstein & Baym 1992, Jones 
1998). In the case that the interaction potential is attractive, 
the coupling always becomes strong at sufficiently small v n 
so that the energy transfer to Kelvin waves must be calcu- 
lated non-perturbatively. In the strong-coupling region, the 
vortex is trapped in the potential for a time t m oc v~ 2 and 
the energy transfer is E a oc v^ 1 . We refer to Jones (1998) 
for details and for a diagram showing a typical vortex orbit 
within the potential. (For weak coupling, realized only at 
values of v n larger than those considered here, the energy 
transfer changes to E a oc v^ 1 ^ and the vortex orbit is little 
perturbed by the potential.) Fig. 1 shows schematically the 
effect of a sequence of such interactions, with a dilute distri- 
bution of point defects, on the vortex path. Each interaction 
produces a lateral displacement and introduces a trapping 
time t m . It is found that, averaged over impact parameter, 
the components of the computed impulse on the vortex are 
approximately equal, 



Fi(£)dt « / F 2 (t)dt 

-oo J — oo 



(1) 



where the axes are as labelled in the Fig. 1. The connection 
between F, defined for interaction with an isolated point 
defect, and the force per unit length f, which is an aver- 
age over macroscopic lengths and times for motion through 
a dense system of point defects, is a complex problem not 
considered here. On microscopic length and time-scales, the 
vortex has fluctuating motion and, as shown in Fig. 1, 
the macroscopically-averaged vortex velocity u = -uu is 
smaller in magnitude than the superfluid velocity v n with 
which the vortex would comove in the absence of interac- 
tion, and has a component perpendicular to it. For general 
values of v n and for a dense system of point defects, it must 
be presumed that as a consequence of equation (1) the force 
f has, in addition to the dissipative component fii , a compo- 
nent fx parallel with k x u, where the circulation is n = kk 
and k is its quantum. In other superfluid problems, the exis- 
tence of this non-dissipative component has been viewed as 
controversial. Thouless, Ao & Niu (1996) have claimed that 
it does not exist for neutral superfluids in a uniform normal 
background (see, however, Hall & Hook 1998, Sonin 1998, 
and the replies therein ), but a system of point defects is not 
of this class. However, in the special case of an infinite linear 
array of point defects, parallel with n, the Green function 
used by Jones (1998) allows one to see that the condition 



F 2 (t')dt' -> 0, 



(2) 



which is necessary for the onset of pinning, can always be 
satisfied at sufficiently small v„. 

The qualitative form inferred for the dissipative force 
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Figure 1. A vortex moves through, and interacts with, a dilute 
system of attractive point defects, shown schematically. In the 
absence of interaction, it would comove with the superfluid veloc- 
ity v n . However weak the attractive potential, the orbit within 
the potential field, at sufficiently small v ni is always character- 
ized by a mean lateral displacement and by trapping for finite 
times in regions close to each origin inside the small rectangu- 
lar boxes. Averaged over macroscopic time and length-scales, a 
vortex typically has velocity vj,, as shown. 



per unit length of vortex, is shown in Fig. 2. Finite 

temperatures and quantum-mechanical tunnelling lead to 
vortex creep and therefore a large and rapidly varying gra- 
dient close to u = 0. However, those properties of fu which 
are important here can be expressed in the simple function 

/||M = Ml-e-/ a )(l-e- 6 /") (3) 

where a, b and f m are constants. We assume that b is of 
the order of v m = fm/pnK, where p n is the density of the 
superfluid continuum outside the nuclei. It is the only nat- 
urally occurring temperature-independent scale of velocity 
in the problem. The precise value of the parameter a is not 
important here except in connection with rotational noise 
generation as mentioned in Section 3; it is very many orders 
of magnitude smaller than b. 

The linearized equations for the macroscopic hydrody- 
namics of an incompressible superfluid rotating with angular 



velocity il are: 

p n K x (v L - v n ) + /||(w)u + f±(u)k x u = 0, (4) 

—7?j- + 2fl x vl = — —VP, (5) 

at p n 

where P is the pressure. Comparison with the equations 
given by Sonin (1987; equations 4.36 - 4.38) shows that the 
Tkachenko force, derived from the vortex lattice structure, 
and the force dependent on vortex curvature have been ex- 
cluded from equation (4). Both are quite negligible in the 
present context. These equations are macroscopic in that 
equation (5) contains the vortex density 2fi//t, although 
equation (4), with no Tkachenko force, is valid for an isolated 
vortex. At the low vortex densities (10 3-5 cm -2 ) of rotating 
neutron stars, v n does not change rapidly with time and, for 
the purpose of understanding vortex unpinning and repin- 
ning at a qualitative level, can be treated as a constant in 
equation (4). Remembering that u = v — vl, and that all 
the velocity components with which we are concerned are 
perpendicular to k, the components of equation (4) parallel 
with u and k x u give an equation for u, 

(v-v n f=(u-±-) 2 + (M\ (6) 

The roots of this equation are determined by the u- 
dependence of the functions fu and f±. For consistency with 
equation (2), the latter must satisfy the limit f±(u) — > as 
u — > 0. In general, there can be either one or three roots 
in each of the sectors | v — v n |< v m and j v — v n |> v m . 
A simple example is the instance in which the solid is sta- 
tionary in the rotating coordinates so that u = — V£, with 
force components f± — 0, and fu given by equation (3). The 
solutions are shown in Fig. 2 for a set of four values of v n . 
In case (1), the single root corresponds with an unpinned 
vortex moving with a velocity vl ~ v n . If v n were decreased 
to case (2), there would be three roots including a pinned 
solution at very small u. Further reduction to case (3) leads 
to a critical value u = u c (u c = 0.44u m for b — 0.3v m ) 
and to a transition to the pinned solution. At the smaller 
v n of case (4), only the pinned solution exists. If equation 
(6) has only one root, there is, of course, no discontinuity 
and u varies continuously with v n . Obviously, this discussion 
is based on our definition and use of macroscopic averaging, 
but gives some idea of the unpinning and repinning processes 
for an isolated vortex. The very low vortex densities in neu- 
tron matter suggest that these, rather than multiple-vortex 
processes involving vortex-vortex interactions, are the sig- 
nificant consideration. 

2.2 The hydrodynamic equations 

Although the non-dissipative force component f± is ex- 
pected to be finite for general values of u, any attempt to 
give it an explicit form would be much less well-based than 
equation (3). Thus the assumption f± = adopted here 
is the only practical basis for obtaining numerical solutions 
of the hydrodynamic equations. Resolved into components 
parallel with the unit vectors f and in cylindrical polar 
coordinates, equations (4) and (5) are, 

/ll,. — p n K (v L 9 - v n e) = 0, (7) 
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Figure 2. The force per unit length fu (in units of f m ) is shown 
versus u (in units of v m ). The function is that of equation (3) 
except that the value of the parameter a has been greatly exag- 
gerated for ease of representation. For the same reason, the second 
parameter is b = 0.3v m , a value intermediate between those con- 
sidered in Figs. 3 and 4. Possible solutions of equation (6), for the 
case v = fj_ = 0, are given by the intersection of this function 
with the circles of radius v n /v m . For v n > v m (case 1), the vortex 
approximately comoves with the superfluid. With decreasing v„ 
(cases 2 and 3), the critical value u c is approached, leading to 
vortex repinning (case 4) 



f\\S + PnK {VLr ~ V nr ) = 0, 



= 2Qv L6 V r P, 

Pn 



V„8 = —2Q.VLr 



— V e P. 

Pn 



(8) 



(9) 



(10) 



With regard to boundary conditions, the only practical way 
of obtaining numerical solutions is to assume an infinite sys- 
tem. Thus the boundary condition that v n should have no 
perpendicular component on the spherical surface of the su- 
perfluid at p n d is neglected. Solutions are therefore valid only 
for wavenumbers q ~ 10 -4 cm -1 large enough to be consis- 
tent with the infinite system assumption. Equations (7)-(10) 
also neglect fluctuations in pressure and gravitational poten- 
tial, a reasonable approximation at the wavenumbers con- 
sidered here. Any pressure fluctuation caused by the local 



shear displacement is negligibly small and it is therefore ap- 
propriate to adopt the initial conditions V r P = VeP = 0. It 
will be convenient to eliminate vl from equations (7)- (10) 
in favour of u. We also neglect any small differences between 
the local k and the z axis. The solid velocity v = v sin tot 
is resolved into components parallel and perpendicular to 
the z axis, only the latter, vi being present in equations 
(7)-(10). The velocity variables are then u r ,ue,v n r and v n e, 
with initial conditions of v nr {0) = and a specified value of 
v„e(0). The initial values, u r ,e{0), obtained from equations 
(7) and (8) are for a pinned-vortex solution which is always 
assumed as the initial condition. Algebraic solution of equa- 
tions (7) - (10) is possible through linearization, which is 
valid for shear waves of such small amplitude that they can- 
not cause vortex unpinning. It gives a dissipation rate per 
unit length of vortex of the order of 



(p n KV j_) 



du 



For any physically reasonable values of the parameters con- 
cerned, this is many orders of magnitude too small to be of 
interest in the present context. 

The rate of energy dissipation per unit length of vortex 
is f • u, which satisfies 



f • u = f • v — f ■ V£. 

Because f + fju 
rewritten as 

f • u = f • v — f • vy. 



and fjf • (vl — v„ 



0, this can be 



(11) 



This merely states that the dissipation rate is the sum of the 
energy loss rates for the shear mode (f ■ v) and for the su- 
perfluid. For the chosen initial conditions, equations (7)- (10) 
have been solved numerically, with attention to the instabil- 
ities discussed in Section 2.1, for several complete cycles of 
the shear mode in order to obtain values of its time-averaged 
energy loss rate. These rates have been calculated for the in- 
terval 10 4 < v a ± < 10 6 cm s _1 , where v x is the component 
of the displacement velocity v perpendicular to k, and for a 
distribution of values of its r, 9 components. This allows us 
to obtain, for a given value of v a and for each of the shear 
mode polarization states, the energy loss rate averaged over 
an isotropic distribution of its wave vector q. (The multiple 
reflection of shear waves makes this average the most ap- 
propriate way of presenting the data.) The energy loss-rates 
per unit length of vortex, divided by p„K and averaged over 
both polarization states, are defined as D av and shown in 
Fig. 3 as functions of v for pairs of values of fm/p n K and 
of b. They are quite insensitive to ui, SI and a. Therefore, 
we have assumed uj = 10 4 rad s _1 , Q, = 10 rad s _1 and 
a — 10~ 4 cm s _1 throughout. Except for velocities close to 
the unpinning threshold, loss-rates do not depend strongly 
on the initial pinned-vortex velocity. This is demonstrated 
by the broken curve in Fig. 3 which is for « r (0) = 0, an ini- 
tial state of an exactly corotating superfluid, v nr = v n e ~ 0. 
Away from the threshold, it is almost indistinguishable from 
the curve for u r (0) = a and so we have assumed this lat- 
ter initial condition for all other curves shown. In general, 
energy loss-rates become large at values of v sufficient to 
produce vortex unpinning and increase to asymptotic values 
pnKD av « / m 6. The mode decay time r at a given ampli- 
tude v can be found from its energy density IpjvWo and the 
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Figure 3. The shear-mode energy loss-rate per unit length of 
vortex, averaged over polarization and for an isotropic distribu- 
tion of wavevectors, is shown as a function of velocity amplitude 
v for given v m and b. The broken curve, for b = v m = 10 cm 
s~ x , is for the initial condition v nr = v n g = 0. 



Figure 4. The vertical axis gives the right-hand side of equation 
(12), which is proportional to Qt averaged over polarization and 
for an isotropic distribution of wavevectors. The parameter sets 
are those of Fig. 3. Damping is strongest at values of v in the 
vicinity of v m . 



loss-rate per unit volume, equal to the product of p n ^D av 
and the vortex number density 2Q/k. The function 

_^£_ = 4Q T £lL (12) 

is shown in Fig. 4. Assuming the lattice parameters calcu- 
lated by Negele & Vautherin (1973), the density ratio p n /pN 
appearing in our, perhaps naive, model of the shear wave 
varies from zero at the threshold p n d to approximately 8 in 
the interval 10 13 < p < 10 14 g cm -3 which includes most of 
the neutron-drip volume. The product ilr is less than unity 
for a substantial interval of shear-wave amplitude. 



3 DISCUSSION AND CONCLUSIONS 

The motivation for expressing the force parameter f m as in 
Figs. 3 and 4 is that v m = fm/pnK is an upper limit for the 
superfluid velocity v n g relative to coordinates rotating with 
angular velocity Q,. Thus values of the order of 10 4 — 10 6 cm 



s _1 must exist in some intervals of matter density if angular 
momentum transfer from crust superfluid is the origin of, or 
a substantial contributor to, pulsar glitches. Much smaller 
values may exist in other regions, and the limit v m — > 
must be approached near p nt j and p c . In Section 2.1, it was 
argued that b should be of the order of v m , but in Figs. 
3 and 4, results are also given for values an order of mag- 
nitude smaller. From Fig. 3, we can see that the asymp- 
totic value D av ~ v m b gives a very approximate estimate, 
2Q.p n Vmb6{vo — Vm), for the energy-loss rate per unit vol- 
ume, where 8 is here the unit step function. The order of 
magnitude of the energy decay time is then 

T _ PnVq , 13 v 

Ap n Q.V m 6(v - V m ) ' 

for b = v m . It is a function only of v /v m , so that even a 
small-amplitude shear wave is strongly damped in regions of 
weak vortex-pinning. It is also temperature-independent, be- 
low the superfluid critical temperature, because it depends 
principally on fn at values of u many orders of magnitude 
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greater than a. The conventionally defined shear viscosity r\ 
for the solid (see Landau & Lifshitz 1970, p. 155) is related 
with the decay time by r\ = p/ruj 2 , where p, is the shear 
modulus. 

This process has some similarity with stress-induced 
motion of dislocations, considered briefly by Blaes et al. 
Both involve dissipative motion of a topological defect, ei- 
ther in the superfiuid order parameter or in the lattice, but 
it is not possible to obtain even an order of magnitude esti- 
mate of the importance of the latter process. The force per 
unit length on a line dislocation is of the order of bpv ±/c s , 
where b is the Burgers vector (see Landau & Lifshitz 1970, 
p. 133); it is typically some orders of magnitude smaller than 
pn^Dav The dislocation density is not known, but could be 
many orders of magnitude greater than the vortex density, 
20/k; « 10 3 " 5 cm -2 . But the time-averaged velocity of an 
individual dislocation under stress is probably very many 
orders of magnitude smaller than v m . 

The coupling with longitudinal sound referred to in Sec- 
tion 1 introduces additional damping. The most significant 
source is the bulk viscosity produced by the strangeness- 
changing four-baryon weak interaction. The amplitude de- 
cay time is 



at matter density p and angular frequency u>, where ce(0) is 
the longitudinal sound velocity in the zero- frequency limit. 
The bulk viscosity coefficient £ is a function both of temper- 
ature and angular frequency. It also depends on the equation 
of state at core densities. For recent evaluations, based on 
differing assumptions, we refer to Jones (2001), Haensel et 
al (2002) and Lindblom & Owen (2002). For u> = 10 4 rad 
s _1 , a temperature of 10 9 K and ( = 3x 10 28 g cm _1 s~ 1 (see 
Jones 2001; equation 45), the decay time is r ni = 3 x 10~ 2 s, 
which would increase by several orders of magnitude if the 
temperature were to decrease to 10 8 K. The evaluation of 
Lindblom & Owen gives a larger £ and thus an even shorter 
decay time. 

Fig. 3 of Blaes et al shows the computed Alfven wave lu- 
minosity as a function of time following generation of shear 
waves by a starquake, assuming a magnetic flux density 
equal to their reference value of 10 11 G. Energy begins to 
enter the magnetosphere about 10~ 3 s after the event and 
the luminosity persists for times just exceeding 6 x 10 _1 s. 
Anomalous X-ray pulsars (AXP) and soft gamma repeaters 
(SGR) are associated with long rotation periods, f2 « 1 rad 
s _1 and, in the magnetar model, with surface magnetic fields 
exceeding 10 14 G. These very large fields indicate, assuming 
the B-dependence of r s given by Blaes et al and referred 
to in Section. 1, that superfiuid friction damping times are 
probably rather longer than the time r s for energy trans- 
fer by coupling of the shear modes with magnetospheric 
Alfven waves. In magnetars, superfiuid friction is unlikely 
to limit the efficiency of this process. A similar conclusion, 
though depending on temperature, holds for the bulk vis- 
cosity damping. 

Duncan (1998) has examined the excitation and de- 
tectability of global seismic modes, particularly low-order 
toroidal modes, by starquake events. He assumes that these 
modes are damped through their interaction with fault lines 
in the solid and the consequent generation of high-frequency 



shear waves. But the superfiuid friction calculated in the 
present paper damps such modes, based on shear defor- 
mations, directly. The damping times found here are much 
shorter than those considered by Duncan and would rule out 
the existence of other than small mode-amplitudes. 

The presence of Q in equations (12) and (13) repre- 
sents no more than the mean vortex density of the super- 
fluid. Therefore, superfiuid friction would produce extremely 
strong damping in more rapidly rotating neutron stars, such 
as the typical radio pulsar with Q, w 10 1-2 rad s _1 . This 
may also be of relevance to the AXP and SGR sources, 
of which all examples so far observed have long periods. 
Energy transfer through coupling between shear waves and 
magnetospheric Alfven waves, as in the magnetar model of 
Thompson & Duncan (1995), is unlikely to be efficient at 
short periods owing to the very strong damping. 

Superfiuid friction is also likely to be relevant to other 
current problems such as the possibility of vortex unpin- 
ning in precessing neutron stars (Link & Cutler 2002), and 
the damping of r-modes in rapidly rotating neutron stars 
(Lindblom, Owen & Morsink 1998; Andersson, Kokkotas & 
Schutz 1999). Many papers on the latter problem have as- 
sumed that the solid crust does not participate in the dis- 
placement amplitude of the mode, is stationary in rotating 
coordinates, and therefore acts only as a boundary condi- 
tion on a thin viscous or turbulent dissipative boundary 
layer. But more recent work (Levin & Ushomirsky 2001; 
Yoshida & Lee 2001) has recognized that the mode ampli- 
tude must penetrate the crust significantly at high values 
of fl. (The radial component of the r-mode displacement, 
though small, must penetrate the crust. Also, an element 
of solid crust moving with velocity v in the rotating coor- 
dinates experiences a Coriolis force which, for fi « 10 2 rad 
s~\ is of the same order of magnitude as the shear stress 
integrated over its surface.) A mode with amplitude large 
enough for its velocity of displacement to be of the order of 
v m will be subject to strong dissipation in the crust, the en- 
ergy loss-rate there being of the order of 2Qp n v m b6(v — v m ) 
per unit volume, as for shear waves. Unlike kinetic theory or 
weak-interaction mechanisms of viscous dissipation (Jones 
2001; Haensel et al 2002; Lindblom & Owen 2002), this is 
not cut off by exponential factors involving superfiuid energy 
gaps but is approximately temperature-independent below 
the critical temperatures. However,the detailed calculation 
required to obtain the r-mode lifetime is beyond the scope 
of the present paper. 

The second term on the right-hand side of equation (11) 
is the energy- loss rate for the neutron superfiuid. For all the 
cases in Fig. 3, except for that with the initial condition 
u r (0) = 0, the superfiuid loses kinetic energy in the rotating 
frame of reference and is gradually brought into corotation 
with the charged components of the star. The initial rate is 
of the same order as for the shear mode. 

Superfiuid spin-down in the presence of shear waves is 
governed by the form of /ii(«) at u < u m , the position u m 
of its maximum being u m ~ Vab for the specific form given 
by equation (3). In the absence of shear waves, the solid 
matter velocity is v — and steady-state superfiuid spin- 
down is possible, in principle, if the required vortex-creep 
velocities, vl = — rtl/2Q, at radius r, are smaller than u m . 
The actual form of /y in this region is unknown, because it is 
dependent on thermal activation and quantum-mechanical 
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tunnelling, but it would not be unreasonable to assume that 
it allows this condition to be satisfied given the very small 
values, vl < 1CT 5 cm s _1 required. Under this assumption, 
some perturbing movement of the solid is needed to induce 
vortex unpinning. Its magnitude depends on the actual form 
of near its maximum at u m , but it is likely that quite 

small values would be sufficient to satisfy the condition | 
v — v„ > v m in one or more, possibly small, regions of 
superfluid. Following our discussion at the end of Section 
2.1, we can see that vortices in these regions would remain 
unpinned with rapid superfluid spin-down, even if the initial 
disturbance were to decay to negligible amplitudes, until the 
critical value u — u c is reached. Thus a significant fraction, 
1 — y '(uc/vm) 2 + {f\\(u c )/fm) 2 , of the available superfluid 
angular momentum in these regions is transferred to the 
charged components of the star in discontinuous glitch or 
noise events. 

It is worth comparing the crust movements considered 
here with those assumed in the glitch model of Ruderman, 
Zhu & Chen (1998). The poloidal velocities of the latter 
model are not necessarily large enough to generate a signif- 
icant shear-wave energy density; it is possible that they are 
no more than episodes of plastic flow occurring in times 
smaller than the observed upper limit (w 10 2 s) for pul- 
sar spin-up, thus with velocities many orders of magnitude 
smaller than the shear- wave velocity c s . In this case, super- 
fluid is spun up or down by the inward or outward radial 
displacement of the solid with pinned vortices (see Jones 
2002). Vortices are not necessarily unpinned although this 
could occur in any region where the superfluid velocity v n 
reaches v m . But crust movement with poloidal velocities of 
the order of c s is not excluded, in which case glitches would 
be initiated by the (v 7^ 0) velocity perturbations considered 
here. It is possible that both kinds of movement contribute 
to the observed glitch population. 

This would not be inconsistent with the value of the 
glitch activity parameter defined by Lyne et al (2000) , which 
is based on the observation that for all except very young or 
old pulsars, the summed glitch amplitudes &Q. g i divided by 
the total observation time is T _1 ^\ 8Q. g i w —1.7 x 10 _2 fi. 
For this group of pulsars, angular momentum of the same 
order as that of the total crust superfluid is transferred by 
glitches rather than by the steady-state superfluid spin-down 
described earlier in this paragraph. Thus all glitches in radio 
pulsars may be initiated by crust movements, with poloidal 
velocities of the order of c s or much smaller. Such movements 
occur less frequently, if at all, in very young or old neutron 
stars for which steady-state superfluid spin-down is more 
usual. These modes of glitch formation may well extend to 
very small events, so giving a model for the observed timing 
noise in the residual phase and its time derivatives. 
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